Stability of off-axis motion for intense particle beams in periodically focusing channels 
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A general equation for the centroid motion of free, continuous, intense beams propagating off-axis 
in solenoidal periodic focusing fields is derived. The centroid equation is found to be independent 
of the specific beam distribution and may exhibit unstable solutions. A new Vlasov equilibrium 
for off-axis beam propagation is also obtained. Properties of the equilibrium and the relevance of 
centroid motion to beam confinement are discussed. 
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A fundamental understanding of the kinetic equi- 
librium and stability properties of high-current, low- 
cmittance beams in periodically focusing systems is cru- 
cial for the development of a wide range of advanced 
particle accelerator and coherent radiation source appli- 
cations. For a long time, the Kapchinskij-Vladimirskij 
(KV) distribution jj was the only Vlasov equilibrium 
distribution known for the propagation of periodically 
focused intense particle beams. Equilibrium and stabil- 
ity analysis based on the KV beam have been critical to 
the development and understanding of the physics of in- 
tense beams @,HSI!H0|ll!|l3. M ° r6 recently ' [t 
has been shown that the KV distribution can be general- 
ized to allow for rigid beam rotation with respect to the 
Larmor frame in periodic solenoidal focusing fields [Tl| . 
Studies indicate that rotation may have an important 
role in particle beam stability [l2| . 

In the derivation of these Vlasov equilibria it is always 
assumed that the beam is perfectly aligned with the sym- 
metry axis of the focusing field [J, H3 ■ Actually, this 
simplifying assumption is generally used in the analy- 
sis of intense beams because the axis is an equilib- 
rium for the beam centroid, and the equilibrium is stable 
if smooth-beam approximations are employed where the 
periodic fluctuations of the focusing field are averaged out 
|l3|. In some cases, however, we may expect the onset of 
parametric resonances involving the centroid motion and 
the focusing field oscillations, which would destabilize the 
centroid motion and heavily affect the overall beam dy- 
namics. In such conditions the averaging procedure is 
no longer valid and a detailed description of the centroid 
dynamics becomes mandatory. 

In this paper, we derive from a kinetic Vlasov-Maxwell 
description a general equation for the centroid motion 
of free, continuous, intense beams propagating off-axis 
in solenoidal periodic focusing fields. It is shown that 
the centroid obeys a Mathieu type equation. The equa- 
tion is independent of the specific beam distribution and 



becomes unstable whenever the oscillatory frequency of 
the centroid, which is related to the rms focusing field 
strength per lattice, is commensurable with the focusing 
field periodicity itself. In the particular case of a uni- 
form beam density around the beam centroid, we show 
that there exists a self-consistent Vlasov equilibrium dis- 
tribution for the beam dynamics. The beam envelope 
that determines de outer radius of the equilibrium beam 
around the centroid is shown to obey the familiar en- 
velope equation [jj 0, H3. Il5l | , being independent of the 
centroid motion. An example of the Vlasov equilibrium 
is discussed in detail to show the possibility of finding 
beam solutions for which the extensively studied enve- 
lope equation |3J, y, |8|, |9|, liil UM * s stable, whereas the 
centroid motion is unstable, revealing the importance of 
the centroid motion to overall beam confinement proper- 
ties. 

We consider a free, continuous charged-particle beam 
propagating with average axial velocity /3(,ce z through a 
periodic solenoidal focusing magnetic field described by 



B(r, s) = B z {s)e z - -B' z (s)e r , 



(1) 



where r = xe x + ye y , r = (x 2 + y 2 ) 1 ^ 2 is the radial 
distance from the field symmetry axis, s — z — flbct is 
the axial coordinate, B z (s + S) = B z (s) is the magnetic 
field on the axis, the prime denotes derivative with re- 
spect to s, c is the speed of light in vacuo, and S is the 
periodicity length of the magnetic focusing field. Since 
we are dealing with solenoidal focusing, it is convenient 
to work in the Larmor frame of reference pj , which ro- 
tates with respect to the laboratory frame with angu- 
lar velocity fiz(s) = qB z (s)/2ji,mc 1 where q, m and 
76 = (1 — P 2 )^ 1 / 2 are, respectively, the charge, mass 
and relativistic factor of the beam particles. The Larmor 
frame is used throughout the paper, such that e x and ej, 
are assumed to be versors along the rotating axes. In 
the paraxial approximation, the beam distribution func- 
tion fb(r, v, s) evolves according to the Vlasov-Maxwell 
system Q 
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— + v • V/fc + {-K z r 
os 

vV = - 



- vvo • v v / 6 = 0, (2) 
(2nK/N b ) n b {r, s), (3) 
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n b = fbdv, 



(4) 



where n b (r,s) is the beam density profile, n z (s) = 
q 2 B 2 (s)/^l(3 2 m 2 c 4 is the focusing field parameter, 
K = 2q 2 N b / r y b (3 2 mc 2 is the beam perveance, N b = 
J f b dvdv =const. is the number of particles per unit 
axial length, and v = r' 17]. In Eqs. ©-Q, i> 
is a normalized potential that incorporates both self- 
electric and self-magnetic fields, E s and B s . It is 
related to the self-scalar and self-vector potentials by 
4> s = {3 b 1 A" z = ^■mP 2 c 2 tp(r,s)/q, where A s (r, s) = 
A s z (r,s)e z , E s (r,s) = -V</> s (r, s), and B s (r,s) = V x 
A-(r.-:, 

Our first task here is to determine the evolution of the 
beam centroid located at 



r(s) = N b 1 J r/ h (r, v, s)drdv. 



(5) 



In order to do that one multiplies Eq. J2J) by r and inte- 
grates over phase-space to get 



r =v, 



(6) 



where v = N b 1 J vfdrdv. If one now multiplies Eq. © 
by v and integrates over phase-space, one obtains 



v' = — K z r — Vip, 



(7) 



where Vip = N^ 1 f \7tpfdrdv is obtained by integration 
by parts of the V v - term in velocity space. Using Eqs. 
© and (@J we can rewrite as 



W> = (2-kK)- 1 J VVA7 2 V>dr. 



(8) 



Then we note that the integrand of Eq. (JHJ can be cast 
into the more suitable form 



VV>V 2 V> = V • [WW - I(V^) 2 / 2 ] 



(9) 



where the unit dyadic I reads I = e x e x + e y e y . Now, 
employing Gauss theorem we obtain 

W = {2ttK)- 1 j> e„ • [VipVif) - I(VV0 2 /2] dA = 0, 

(10) 

because Vi/> — > as r — > oo for beams in free space. 
In Eq. lfTU|) . dA and e„ are, respectively, the boundary 
differential element and the unit vector normal to the 
boundary of integration located at r — > oo. In fact, the 
result Vip = is expected based on the action-reaction 
law; since — Vip corresponds to the self-force exerted on 
the beam particles by themselves, its average throughout 
the beam distribution has to vanish due to the pairwise 
structure of the interparticle electromagnetic interaction. 
Using Eqs. Q and (fTU|l in Eq. J(jJ), we finally obtain the 
centroid equation of motion 



Let us stress that we have not made any assumption 
on the particular form of the beam distribution func- 
tion so far. Thus, the centroid equation above is al- 
ways valid as long as the beam evolves according to the 
Vlasov-Maxwell system, Eqs. (0)-(@J). In the laboratory 
frame, combined to the oscillatory motion described by 
Eq. (|ll|l the centroid also rotates with angular veloc- 
ity f2i(s) around the center r = 0. Taking into ac- 
count that k z (s) is periodic, Eq. (|llfl is of the Math- 
ieu type which is known to present unstable solutions 
related to parametric resonances in the r motion. If we 
conveniently write the average of k z (s) over one lattice 
period as (1/S) n z {s)ds = <j 2 /S 2 , where cr n is a dimen- 
sionless parameter proportional to the rms focusing field 
measuring the vacuum phase advance in the small field, 
smooth-beam approximation, the instabilities in the cen- 
troid motion are expected when one approaches <7o ~ 7w; 
this condition corresponds to parametric resonances be- 
tween the oscillation periodicity of f in the average (rms) 
focusing field and the periodicity of the focusing field it- 
self. Depending on the exact profile of k z (s) the size of 
the unstable regions surrounding Co ~ nn vary signifi- 
cantly. If the aim is beam confinement, these regions are 
to be avoided. 

It is worth mentioning that although Eq. I|ll(l is 
strictly valid for free beams only, it is expected to pro- 
vide a good description of the centroid motion in bounded 
systems if the beam is nearly symmetric and is not ex- 
cessively displaced from a pipe center located at r = 0. 
The reason is because in this case V-0 = ±e„|V-0| at the 
pipe walls, where e n is now the unit vector normal to 
the wall, and the surface integral in Eq. QlUfl. performed 
along the boundary, still vanishes since |V^| is approx- 
imately constant there. Note also that the presence of 
a pipe would generally not suppress the centroid insta- 
bilities discussed in connection with Eq. (jllfl : in fact, it 
would even enhance it because the image charges induced 
are of opposite sign, attracting the beam to the wall. 

Our next task is to show that we can construct a Vlasov 
equilibrium for off-axis beam transport. In particular, we 
assume a beam with a uniform radial density distributed 
around a center located at r (s) = x (s)e x +j/ (s)e y , i.e., 



n b (r,s) = 



N b /nr 2 b (s), 
0, 



r s < r b (s), 
r s > r b (s), 



(12) 



where r b (s) is the equilibrium beam envelope and = 
r r Q . A schematic of the beam distribution of Eq. i|12|) 
and corresponding vectors is shown in Fig. 1. For such 
beam we can easily recognize r D (s) as being the centroid 
coordinate. According to what was shown previously, its 
evolution must then obey 



r o " + K z (s)r o = 0. 



(13) 



Using the prescribed n b (r,s) in Eq. © we find for the 
normalized self-potential 



+ K z ( S )f = 0. 



(11) 



4,(v, S ) = -Kr 2 /2r 2 {s) 



(14) 
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in the beam interior (r$ < r&). Therefore, a single par- 
ticle of the beam located at r(s) subjected to the exter- 
nal focusing field force — n z (s)r and the self-field force 
— V?/>(r, s) will evolve according to 

v" + K z v-{K/rl) n = 0. (15) 

If we now subtract Eq. from Eq. (|15|) we obtain 

y s " + n z v s -{K/r 2 b )v s = 0, (16) 

which describes the motion of the beam particle with 
respect to the center r Q . Equation l|16f) can be solved 
with known techniques of physics of beams 0, Q . Con- 
sidering the motion along the x-axis, we write x$ — 
A xs w(s) sin[/ s ((s)ds+( xs o] with A xs and ( xsQ constants. 
Substituting this expression into Eq. ^ we obtain 
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k(s)w 



w 



(17) 



C(s) = w 2 (s), where k(s) = k z (s) — K/r 2 {s), and the 
constant of motion A xs can be expressed in the form 



A,,, 



(xs/w) 2 + (wx' s - w'xsY 



(18) 



Performing an equivalent calculation for the motion along 
the y-axis, one shows that A yg given by 



A t s = {ys/wf + {wy' s - w'ysY 



(19) 



is also a constant of motion. From Eq. (|16[) one sees 
that all the forces are central with respect to the centroid 
r Q . Thus, one readily demonstrates that the canonical 
angular momentum Pqs given by 



Pes = xsy' s - ysx' s 



(20) 



is a constant of motion as well. Because A 2 S , A 2 6 , and 
Pqs are exact single-particle constants of motion, a pos- 
sible choice of Vlasov equilibrium distribution function 
is 



if Q (r,v, S ) 



— — ' 

7T Z €T 



A,. 



A - 

ys 



2w h Pe<5 - (l - u> 2 ) e T 



(21) 



where dfF® /ds = 0, ct =const is an effective emittance, 
and the rotation parameter LUb =const is in the range 
— 1 < Ub < 1 for radially confined beams. Using f^ in 
Eq. (0J , it is readily shown that the uniform density pro- 
file centered at r Q of Eq. (|12|) is consistently obtained, 

1/2 

provided n,(s) = e T ' w(s). Hence, rt(s) obeys the famil- 
iar envelope equation 



K 



Kz(s)r b j = 0. 



n 



(22) 



Performing the appropriate averages over the equilibrium 
distribution, Eq. I|21|) , we can show that the beam rigidly 
rotates around its centroid r Q (s) with angular velocity 



FIG. 1: Schematic of the beam distribution of Eq. 1121 and 
corresponding vectors. C corresponds to the centroid posi- 
tion. 



^bs(s) — ujb^TPbc/r 2 {s). Also, the rms emittance calcu- 
lated with respect to the centroid position is given by 



<x 2 ) <x' 2 ) - (x xf 



1/2 



(23) 



where \ — xs,ys, and the brackets indicate averages over 
the beam distribution. One thus sees that a Vlasov equi- 
librium distribution can be formed for which the beam 
envelope obeys Eq. lll'L'1) with constant emittance even 
when the centroid moves off-axis, r Q ^ 0, following the 
dynamics dictated by Eq. (JT3J. We refer to this equilib- 
rium as a periodically focused off-axis Vlasov equilibrium. 
Let us call attention to the interesting fact that the cen- 
troid motion and the envelope dynamics are uncoupled in 
this case. In other words, centroid dynamics does not af- 
fect the known stability results for the envelope dynamics 
LH S El EH and is not affected by the latter as well. 
One should keep in mind that for good beam confinement 
both centroid and envelope have to be stable. 

We now illustrate our results with an example of peri- 
odically focused off-axis Vlasov equilibrium. We consider 
a particular set of parameters for which the envelope 
equation ll'L'l) is known to be stable, whereas the cen- 
troid motion of Eq. I|13|) was found to be unstable. We 
investigate beam transport with the aid of self-consistent 
numerical simulations, where a large number Nb — 8000 
of macroparticles interact via pairwise electromagnetic 
interactions 9]. In the simulation we used SK/er = 5.0 
and S 2 k z (s) = cr 2 [l + cos(27rs/S')], with a = 155°, over 
20 lattice periods. The macroparticles were launched at 
,s = according to the equilibrium distribution, Eq. I|21|) . 
with ujb = 0, r Q = = r Q ', and corresponding to the 
matched solution with rb(s + S) = rb(s) of the enve- 
lope equation l|22|) . The finite number of macroparticles 
in the initial condition acts as a seed for any possible 
instability to develop. Simulation results are presented 
in Fig. 2. The evolution of the centroid displacement 
r a = |r Q | calculated from the macroparticles positions 
r as r Q =< r >, where the brackets indicate average 
over macroparticles, is shown in Fig. 2(a) (circles). It 
reveals that the centroid motion develops the typical ex- 
ponential growth of unstable dynamics that agrees with 
the fact that the set of parameters considered leads to 
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FIG. 2: Multiparticle self-consistent simulation results, (a) 
The centroid motion; (b) rms emittance; and (c) the envelope 
dynamics. Centroid displacement and envelope are normal- 
ized to (Se T ) 1/2 . 



an unstable solution of Eq. l|13f) . The solid line corre- 
sponds to the solution obtained from Eq. i|13[l . Despite 
the centroid instability, the beam equilibrium distribu- 
tion is preserved as verified in Fig. 2(b) that shows that 
rms emittance is well conserved as the beam evolves. 
RMS emittance is calculated according to Eq. i|23[l . 
considering averages over macroparticles. Finally, Fig. 
2(c) compares the envelope obtained from the envelope 
equation (|22|l with that obtained from the simulation, 
r b = [2 < (r - r Q ) 2 >] 1/2 , for the last 3 periods of the fo- 
cusing channel. The perfect agreement proves once more 
the preservation of the equilibrium distribution. More- 
over, we see that in spite of the unstable centroid the 
envelope is stable, as predicted. 



To conclude, based on kinetic grounds we have derived 
a general equation for the centroid motion of free, con- 
tinuous, intense beams propagating off- axis in solenoidal 
periodic focusing fields. It was shown that the centroid 
equation is independent of the specific beam distribu- 
tion and may exhibit unstable solutions. In the partic- 
ular case of a uniform beam density around the beam 
centroid, we have shown the existence of a periodically 
focused off-axis Vlasov equilibrium distribution describ- 
ing a beam that rigidly rotates with a prescribed angular 
velocity around a moving centroid. The beam envelope 
around the centroid was shown to obey the familiar enve- 
lope equation, being independent of the centroid motion. 
An example of periodically focused off-axis Vlasov equi- 
librium was discussed in detail to show the possibility of 
finding beam solutions for which the envelope equation 
is stable, whereas the centroid motion is unstable, re- 
vealing the importance of centroid motion to the overall 
beam confinement properties. 
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